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CLh : ABSTRACT 
PP; 

Q_j Several planets have recently been discovered around stars that are old and 

Q . metal-poor, implying that these planets are also old, formed in the early universe 

^ , together with their hosts. The canonical theory suggests that the conditions 

! for their formation could not have existed at such early epochs. In this paper 

! we argue that the required conditions such as sufficiently high dust-to-gas ratio, 

^ I could in fact have existed in the early universe immediately following the first 

! episode of metal production in Pop. Ill stars, both in metal-enhanced and metal- 

! deficient environments. Metal-rich regions may have existed in multiple isolated 

O i pockets of enriched and weakly-mixed gas close to the massive Pop. Ill stars. 

^ ! Observations of quasars at redshifts z ~ 5, and gamma-ray bursts at z ~ 6, show 

^ ! a very wide spread of metals in absorption from [X/H] ^ —3 to ^ —0.5. This 

I suggests that physical conditions in the metal-abundant clumps could have been 

. ^ ! similar to where protoplanets form today. However, planets could have formed 

^ i even in low-metallicity environments, where formation of stars is expected to 

c3 I proceed due to lower opacity at higher densities. In such cases, the circumstellar 



accretion disks are expected to rotate faster than their high-metallicity analogues. 
This in turn can result in the enhancement of dust particles at the disk periphery, 
where they can coagulate and start forming planetesimals. In conditions with the 
low initial specific angular momentum of the cloud, radiation from the central 
protostar can act as a trigger to drive small-scale instabilities with typical masses 
in the earth to jupiter mass range. Discoveries of planets with low-metallicity 
hosts (e.g. HIP 11952 with [Fe/H] ~ —1.95) shows that planets did indeed form 
in the early universe, which may require modification of our understanding of the 
physical processes that produce them. This work is an attempt to provide one 
such heuristic scenario for the physical basis for their existence. 
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Subject headings: planetary systems: formation - quasars: abundances - cosmol- 
ogy: early universe 



1. Introduction 



A two-year gravitational lensing survey by the Microlensing Observations in Astro- 
physics (MOA) and Optical Gravitational Lensing Experiment (OGLE) groups towards the 
Galactic bulge (Sumi et al. 2011) has found ten events which can be attributed to Jupiter- 
mass planets. These microlensing planets are free-floating, in the sense that no host stars 
have been detected within about 10 AU. Recent est imates show tha t there can be up t o 
100,000 as many such planets as stars in the Galaxy (jSumi et al.ll201ll : IStrigari et al.ll2012l ). 
Three out of the 10 microlensing events from planets have galactic latitudes of 6 ^ —3°, 
and one has b ^ —6°, corresponding to the heights from 0.4 to 0.8 kpc. Therefore these 
planets are likely to belong to the thick disk population, where a considerable fraction (up 
to 30%) of stars is o ccupied by the oldest (Pop. II) Milky Way stars with m ean metallic- 
ity ([Fe/H]) = -0.6 dZoccah et ahlboosi liensbv et al.ll201ll : IffiU et al.lboilH ^. Thus these 
planets could also be old, formed in the early Milky Way Galaxy more than 10 Gyrs ago; 
their free-floating status a result of being ejected from the parent systems by planet-planet 
scattering events during the early stages of planetary systems hfe (e. g.. Ford, Rasio & Yu 
2006). 

In general, a well-known planet distribution function peaking at solar metallicity (Udry 
& Santos 2007) is readily understood in terms of the thermodynamics of the gas from which 
the planets are formed: the smaller the gas metallicity, the weaker is the radiation cooling 
and therefore the higher are the gas temperature and the Jean's mass. In the extreme case 
of primordial gas where the metallicity [Fe/H] = — oo, the only source of cooling is from 
molecular hydrogen lines, and only extremely massive stars may be formed, M > 80 Mq 
(eg. Schneider et al. 2003, Bromm & Larson 2004), which would have evolved rapidly 
and be extinct by now. However, it was suggested (Nakamura & Umemura 2001) that the 
fragmentation of filamentary primordial gas clouds could result in the formation of low-mass 
stars even in the absence of metals, with only molecular hydrogen as the main cooling agent. 
Recent simulations (Clark et al. 2011, Greif et al. 2011) demonstrated that the fragmentation 
process is rather environment-sensitive, such that gravitational instability may result in the 
formation of a group of low-mass (M ~ 0.1 — 10 Mq) protostars with a relatively flat 



^[X/H] = logio (Nx/NH)star ~ logig (Nx/Nr)© and the subscript '0' refers to the Sun. 
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initial mass function (IMF). The discovery of the extremely metal-poor (Z < 7.4 x 10~^)g 
low-mass (M < 0.8 Mq) star (Caffau et al. 2012) seems to confirm this conclusion. The 
formation of planetary fragments in gas with primordial chemical composition is, however, 
problematic, particularly because H2 and HD molecules provide quick thermalization of gas 
with the cosmic microwave background, which limits the fragment mass to M ^ 0.1 Mq 



( IShchekinov fc VasilievI 120061 ). 



The observed "metallicity effect" may be a result of a strong selection effect in that 
metal-poor Pop. II stars in the solar vicinity are at least two orders of magnitude less 
populous than those of solar metallicity (Gilmore et al, 1989; Dehnen & Binney 1998). 
Classical detection methods (radial velocity measurements and transits) limit detection of 
planets to nearby stars and thus most stars with detected planets have metallicities close 
to solar (see, for example. Fig. [T]). Moreover, the metallicity effect is essentially limited 
to giant planets, which are obviously easier to detect. The recent population synthesis 
study (Mordasini et al. 2012) have shown that high metallicity favours the formation of a 
larger number of high-mass planets. The current ratio of giant (M > 0.3 Mj) to small- 
mass (M < 0.3Mj)|3 planets is ~ 3.4, which can indicate that more metal-rich hosts are 
overrepresented (see Fig. [2]). 

In the lower mass domain — Neptune masses — this effect weakens (e.g., Mordasini 
et al. 2012). Indeed, host stars of Neptune-mass planets with precisely measured radial 
velocities (Gumming et al. 2008) show apparently a flat metallicity distribution (Udry et al. 
2006; Sousa et al. 2011): a Gaussian peak at [Fe/H] = —0.5 transforms to a wide plateau 
in the range [Fe/H] = —0.9 to —1.5. No metallicity trend is predicted for even lower-mass 
planets (Mordasini et al. 2012); moreover, recent results from Kepler mission (Buchhave et 
al. 2012) demonstrated that formation of small terrestrial-size planets does not necessarily 
require a metal-enriched environment. 

That planets can form from matter with lower than solar metallicity is evident from the 
fact that about 50% of the giant planet-host stars have sub-solar metallicities (Setiawan et 
al. 2012). A number of giant planets at fairly low [Fe/H] is now known (in 2012 there are 27 
confirmed planetary hosts with [Fe/H] < —0.4, 3 of them with [Fe/H] < —1.0, and a pulsar 
B 1620-26 planet in an [Fe/H] = —1.5 globular cluster M4) and this number is constantly 
increasing; the most recent discovery of two giant planets around an [Fe/H] = — 1.9 F-dwarf 
HIP11952 suggests the age of the system of ~ 12.8 Gyrs (Setiawan et al. 2012). 



is the mass fraction of metals 
^Mj is the Jupiter mass 
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It is worth mentioning that when the first population statistics on exoplanets was pub- 
hshed in 2007 (with very few low-metal hosts known), the relative percentage of metal-poor 
stars ([Fe/H] < —0.6) with planets in solar neighbourhood (~ 4%) was four times higher 
than the relative percentage of Pop. II stars (~ 1%) (as can be inferred from the Fig. 9 in 
Udri & Santos, 2007). According to the latest data (November 2012), this percentage now 
is around 2% (8 out of 472 stars with measured metallicity within 600 pc of the Sun have 
[Fe/H] < —0.6, see Fig. 3). With more discoveries, the relative number of metal-poor hosts 
may asymptotically approach the relative number of Pop. II stars in the solar neighbour- 
hood. Thus, the question of principal importance arises of how to collect sufficient amount 
of metals to form planets, is it possible at all, and what could be the scenarios? 

In this note we argue that the conditions for selective assembling of metals in the dust 
form could have existed even in the very early universe immediately after the formation of 
first massive stars, i.e. at the epoch when the first stellar light illuminated the universe at 
redshifts ^ ^ 10 (see, for ex., Scannapieco et al 2006). It is worth mentioning in this connec- 
tion that very recently Wickramasinghe et al. (2012) presented arguments that primordial 
free-floating (rogue) planets of solid hydrogen may account the whole "missing baryons" in 
the Universe. When recalculated for the Milky Way the number of such primordial planets 
is as numerous as ~ 10^^. In this work we however address rather the question of whether 
planets can form around stars with low metallicity. Moreover, we do not consider possible 
scenarios of planet formation in such conditions, but focus primarily on even a simpler issue: 
what mechanisms can provide accumulation of sufficient amount of metals (mostly in solid 
phase) to make the subsequent formation of planets possible within scenarios similar to those 
acting in the local Universe. 

In Section 2 we bring evidence that in the very beginning of stellar nucleosynthesis in the 
universe, one might find enriched regions with the metallicities close to solar, or even higher, 
where planets can form under conditions analogous to the ones in local univers e , prov ided 



that dust - to-me tal ratio is of ~ 0.3 (as is normally assumed, see iMaiolino et al.l (120041 ) and 



Fukugital (l201ll )). In Section 3 we describe a possibility of enhancement of dust abundance 
in metal-deficient gas either due to the centrifugal force in peripheral regions of circumstellar 
protoplanetary accretion disks, when the specific angular momentum is sufficiently high, or 
due to the action of the radiation pressure from the host protostar in the conditions of small 
specific angular momentum. Section 4 summarizes our results. 

We have used for this discussion the observational data from the Extrasolar Planets En- 
cyclopaedia (http : //www . e xoplanet . en) maint ained by J. Schneid er, from the Exoplanet 



Orbit Database (http: //www .exoplanets. org) (jWright et al.ll201ll ) and from NASA Exo- 



planet Archive (http://exoplanetarchive.ipac.caltech.edu). For the planet hosts with 
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Fig. 1. — The plot shows the distribution of stars with confirmed exoplanets with metal- 
hcity (total 676). It can be seen that most of the planetary hosts are located in the solar 
vicinity, within < 600 pc of the Sun. This may indicate the selection effect in the suggested 
metallicity-planet-hosting dependence. This figure was made using the Extrasolar Planets 



Encyclopaedia at http : //www . exoplanet . eu 
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Fig. 2. — The plot shows the distribution of the exoplanets masses vs host metallicity. It is 
clearly seen that most of the detected planets are giants and the observed "metallicity effect" 
is neatly reproduced by giants, i.e., the planets with mass M > 0.3 Mj (Right). In the small- 
mass domain, M < 0.3 Mj (Left), the mass is nearly uncorrelated with metallicity. This 



figure was made using the Extrasolar Planets Encyclopaedia at http://www.exoplanet.eu 
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Fig. 3. — The plot shows the distribution of planets hosts in the solar vicinity (within 600 
parsec of the Sun) as a function of host metallicity. Out of 472 planetary hosts with measured 
metallicity, 8 have [Fe/H] < -0.6 (~ 2%). 

no metallicity data from the above-mentioned websites, we have used the values from Mal- 
donado et al. (2012). 

2. Metallicity in the early universe 

It became clear from the very first detections of metals at high redshifts 2; ^ 6 (Songaila 
2001) that there are extreme variations of metallicity in the universe. Further observations 
confirmed this conclusion. Schaye et al. (2003) demonstrated that the abundance of carbon 
in the intergalactic medium can vary up to 2-3 orders of magnitude, particularly at lower 
redshifts [z = 1.5 — 2.5), while Simcoe et al. (2006) found many clumps and large clouds with 
metallicity varying between [X/H] = —3 and 0.6 toward quasar HS1700+6416 (zabs = 2.73), 
here X denotes mostly the abundances of carbon, silicon, magnesium ions. Schaye et al. 
(2007) inferred a wide range (from [X/H] = —1 to 0.5) of metallicities in the absorption 
systems of several quasars at redshifts of ^abs = 2 to z^hs = 3.2. They concluded that the 
metallicity tends to be higher for smaller sizes of the absorbing gas clouds, as is predicted 
by the numerical simulations of metal- mixing by Dedikov & Shchekinov (2004). 

From the theoretical point of view three basic models explaining the inhomogeneous 
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metal distribution in the universe have been described: patchy (unpercolated) distribution 
of metal-enriched bubbles around active (starburst) galaxies (Nath & Trentham 1997; Gnedin 
1998; Ferrara et al. 2000); dependence of the metallicity of enriched bubbles on the host 
galaxy (Nath & Trentham 1997; Gnedin 1998; Ferrara et al. 2000; Shchekinov 2002), and 
inefficient (incomplete) mixing (Dedikov & Shchekinov 2004). 

More evidence for highly inhomogeneous metal distribution at high redshifts comes from 
observations of Gamma-Ray Bursts (GRB) and particularly their optical afterglows: optical 
absorption spectra of the so-called GRB damped Ly-a systems - as opposed to QSO damped 
Ly-a systems - and the stellar light from GRB hosts show a high spread of metal abundance 
of about two orders of magnitude, with a weak redshift evolution and a relatively high mean 
[X/H] = —0.7 (Savaglio 2008; Savaglio et al. 2009), suggesting that the metal distribution 
is highly inhomogeneous in the host galaxies. Numerical simulations of mixing of metals in 
the ISM of the Milky Way Galaxy (de Avillez & Mac Low 2002) have shown that it can take 
as long as 120 Myr on length scales of 25 to 500 pc. If this time-scale is also applicable to 
GRB hosts, then homogenization of metals would take a considerable fraction of the Hubble 
time at 2; = 8. The high spread of metallicity in the GRB hosts then simply reflects the 
inhomogeneity of their ISM. Under such circumstances, a wide range of physical conditions 
favourable for the formation of planets around hosts with solar or near-solar metallicity can 
be found in the early universe. 



3. Low-metallicity environment 

Although we have shown that high-metallicity regions favourable for forming planets 
may be common in the early universe, the independence or even possible anti-correlation 
of the frequency of small planets with metallicity (Mordasini et al. 2012; Buchhave et al. 
2012), together with the existence of old metal-poor stars hosting planets, such as for ex. 
HIP 13044 ([Fe/H] ^ -2, Setiawan et al. 2010) or HIP 11952 ([Fe/H] = -1.8, Setiawan et 
al. 2012) suggest that metallicity may not be the governing parameter for the formation of 
planets. It seems therefore worthwhile to consider possible mechanisms which can promote 
formation of planets in a metal-poor gas. 

Core accretion scenario is recognized to account for most of detected planets. However, 
recent discoveries of planets which lie far from the possible parameter space for core accre- 
tion, for instance, planetary mass-metallicity relation, indicate that there may be several 
formation scenarios that lead to systems very different from the usual (see e.g., Mordasini et 
al. 2012). Such planets (for example, HR 8799, BD-f20 2457 or HIP 1952) are difficult to ac- 
count for by the core accretion. Planets found around stars with low metallicities or on very 
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wide orbits pose a challenge to the core-accretion scenario, and a gravitational instability 
model was proposed to better explain their formation (Dodson- Robertson et al. 2009). The 
scattering of these planets to wide orbits as a scenario of their formation (within the core 
accretion framework) was also rebutted in Dodson- Robertson et al. (2009). It seems natural 
to assume that stars with different metallicities form their planets by different mechanisms. 
Such an assumption is consistent with the findings that both giant and Neptune-mass planets 
discovered so far have rather flat metallicity distribution (Santos et al. 2004; Udry & Santos 
2007; Sousa et al. 2008, 2011; Mayor et al. 2011). If the dominant planet formation mech- 
anism were not changing with metallicity, we would rather see a monotonically increasing 
detection rate with [Fe/H] for all metallicity range (Dodson- Robinson 2012). 

Using order of magnitude estimate, we argue that, in the process of star formation 
from metal-underabundant material, immediately after the primeval enrichment in the early 
universe, selective assembling of dust around protostars can operate and enhance the local 
dust abundance to a sufficient level so that physical mechanisms (such as dust coagulation) 
leading to planet formation can start operating. At this stage, we believe that such an 
order of magnitude calculation is more practical than a sophisticated numerical approach 
involving the entire evolutionary path from dust coagulation through planetesimals to planet 
formation - a scenario which is still not fully understood even for our own Solar System. 



3.1. Gravitational collapse versus gas metallicity 

Two key ingredients controlling formation of planets are the presence of dust in a star- 
forming medium and the amount of rotation. It was recently recognized (see, for ex. Toddini 
& Ferrara 2001; Nozawa et al. 2003 and Maiolino et al. 2004) that dust particles are injected 
along wit h metals r i ght fr om the very initial metal-pollution episodes. We therefore assume. 



following iFukugital (1201 if ), the mass abundance of dust to be D ~ 0.3Z, where Z = Zq^ is 



the absolute mass abundance of metals. 

It is obvious that in early galaxies with metal-poor gas, the low amount of metals and 
dust restricts the formation of rocky planets. The maximum number of Earth-size planets 
that can form in a protostellar cloud of mass M is K ~ DM /Me where Me is the Earth 
mass. Defining m = M /Mq, we obtain K ~ 3 x lO^m^ per cloud. For the protostellar cloud 
of mass m = 10 and even at 3 ~ 10"^ and efficiency of planet formation of 1%, we obtain 
K ~ 3. If we assume the star formation efficiency at 10% (i.e. each cloud of m = 10 produces 
one star), it would mean that every star may have 3 planets. Note that when applied to 
the cloud with solar metallicity, this estimate gives up to 300 planets per star. Although it 
was recently recognized that multiplanetary systems are rather a rule than an exception, the 



- 9 - 



most populous exoplanet systems yet detected, HD 10180, has 7 confirmed planets (Lovis et 
al. 2011). It is likely that the population of planets in a given system is usually less than 
10 (more likely 3) because of the instability of planetary orbits on time scales greater than 1 
Myr (Juric 2006; Vedas et al. 2011). The detection of multi-planetary systems can, however, 
be impeded by the high inclination of the planetary orbits to the line of sight (Tremaine & 
Dong, 2012) and/or by the weak and tangled wobbling of the host star produced by distant 
planets. 

It is assumed after Hoyle (1953), that during gravitational contraction a cloud radius 
R can be subjected to a successive (hierarchical) fragmentation, if radiative cooling is suf- 
ficiently efficient to keep temperature growing slower than T oc ra^/^, such that the Jean's 
mass Mj ~ 30T^/^?t,~^/^ Mq decreases with gas density n. More recent numerical simulations 
showed that successive fra gmentation co ntinues while gravitational contraction is isothermal. 



i.e. 7 = d\n.p/d\n.p < 1 ( iLi et al.ll2003l ). In general, since optical depth grows as r oc -R^^, 
a contracting cloud becomes eventually opaque, even though at initial stages it is always 
optically thin. This moment finalizes the isothermal regime and determines the minimum 
stellar mass. 

A more recent development of this approach applied to the problem of Pop. II. 5 stars 
(stars form ed right after Pop. Ill) was developed by Omukai et al. (2005), Schneider et al. 



(2006) and [Schneider et all (120121 ). They considered a parcel of gas density p{t) belonging 
to a contracting fragment whose temperature is determined by heating from gravitational 
compression, and cooling is determined by radiation liberated in inelastic collisions, 

1 k dT d 1 

-P—--L, (1) 



7 — 1 pmii dt dt p 

where mn is the hydrogen mass, p is the mean molecular weight and L is the cooling rate. 
Eq. ([1]) determines the evolutionary track T(p). For a given model of dust particles and 
their mass fraction, and when all possible cooling mechanisms, such as cooling in atomic (CI, 
CII, and 01) and molecular (H2, HD, OH, CO) emission lines, cooling in continuum due to 
thermal emission by dust particles and others are included, the minimum mass of protostellar 
fragments is ~ 0.3 Mq for metallicity [X/H] = —7, and ~ 1 M© for [X/H] = (Omukai et 
al. 2005). Therefore, within these assumptions, the masses of protostellar condensations are 
rather insensitive to gas metallicity from which they form, but can be se nsitive to a dust 



mode l and to the depletion factor of refractory elements on dust surface (jSchneider et al. 



2OI2I ). For simplicity, we restrict ourself with the models described by Omukai et al. (2005) 



and Schneider et al. (2006). 

The cloud density at which fragmentation stops is, however, rather sensitive to metallic- 
ity and can be roughly approximated as n ~ lO^j"^'^ — 10^3"^ cm~^, depending on the mass 
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of a progenitor of the polluting supernova (SN), 195 and 22 M©, respectively (Schneider 
et al. 2006). This can mean that less-enriched protostellar clouds and, as a consequence, 
the circumstellar (protoplanetary) disks when established, must rotate faster by factor of 
ra^/^ if they had equal specific angular momentum jo on average. The specific angular mo- 
mentum of star-forming molecular clouds in the Milky Way is jo ~ 30 pc km/sec, nearly an 
order of magnitude low er than the specific angular momentum in the interstellar medium 
( iBlitz fc Williamslll999l ). which we will assume here as a reference value jo = lO^^A cm^/sec, 
with A as a free parameter. 



3.2. Centrifugal assembling of accreting dust 

Planets seem to form during subsequent evolution of such a protostellar cloud, after 
formation of embryo disk through coagulation, settling of dust onto plane and followed 
fragmentation. The Toomre parameter for the stability of disk (Binney & Tremaine 1987), 

seems to be relatively insensitive to metallicity as both angular velocity Q and surface density 
S scale with the radius as R~^. Q is sensitive to the free parameter A, which can vary in a 
wide range. The sound speed Cs at the time when fragmentation stops and accretion disk 
is to be formed depends on meta llicity, c.s oc ?)~ ^ ^\ an d ca n vary from 0.3 t o 3 km /sec for 



3 = 0.1 — 10 ^ (as inferred from lOmukai et al.l (120051 ) and [Schneider et al.l ( l2006l )). Later 



on, when the circumstellar accretion motion settles on to the disk, the sound speed relaxes 
to the value determined by the heating from central star (protostar) and by the viscosity 
in the disk. From this point of view, and particularly accounting for the fact that A varies 
in a wide range and might be low, one can conclude that the conditions for gravitational 
instability in circumstellar disks today and in the early galaxies are similar. 

At the initial stages when the disk has only just formed, dust particles move tightly 
coupled to gas. This is a fast rotation with velocity v,j) = QR and a slow inward accretion 
with velocity Vr ~ av'^/QR, a < 1, where Vd is the dispersion velocity in gas (see Pringle 
1981). Assuming the turbulence in the disk to be subsonic, let f d < oc Dust 
decouples from gas when the coUisional drag rate, 

Ud ~ oc a , (3) 

rrid 

becomes smaller that the rotation rate VL due to the growth of dust particles; and a are 
the grain mass radius, respectively. The total rate of the dust growth Vg through coagulation 
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and freezing of heavy elements on dust surface can be estimated as 

Ida smuV^nZ 

vAa) = --r ~ , (4) 

' adt ap ' ^ ' 

where ~ 3 g/ cm^ is the density of a dust grain and s ^ 1 is the mean sticking coefficient 
in dust-dust and atom-dust colhsions; tight colhsional couphng between dust and gas is 
exphcitly assumed. After time t = 27iN/Q, corresponding to rotational periods of the 
disk, the grain radius can grow up to 

smnnvdZ 2ttN , , 

a z — , 5 

p \l 

where we explicitly assumed that all dust particles are of equal radius. With the growth of 
dust grains, the collisional drag rate Ud decreases as 

so that after > Ud/^ ^ 1, dust particles decouple from the gas motion. Though 

they continue the frozen rotation, the inward radial motion slow down dramatically, Vdr/vr ~ 
Oiyl/v"^), and dust piles up in a ring. This happens when the dust particles have moved 
inwards by 

27TNVr ^7 fn 

Afl.^.2xio-^. (7) 

where m is the mass of a protostellar cloud in solar units. The gas inflow continues, but 
its faster decoupling {AR/R is small) implies that it is left behind at larger radii. When 
the density at last fragmentation episode n oc 3"'' (6 = 1.5 — 2, Schneider et al. 2006) is 
substituted, it becomes readily seen that after decoupling, dust particles in metal-deficient 
disks accumulate mostly in outer regions of the disk as AR oc 3*"^, i.e. the lower the 
metallicity of the disk the farther in the disk dust accumulates, everything else (e.g. A, 
m) being equal. For the abundance pattern and dust yield corresponding to a M = 22 Mq 
progenitor, we numerically obtain 

AR 0.1m2/33i/3 
^ A^' 

which shows that in metal- deficient circumstellar disks dust accumulates in external regions, 
while in the cases with 3 ~ 1 dust occupies nearly the entire disk. Asymptotically, the 
enhancement factor in the region of accumulation is approximately 77 ~ AR/6R; here SR is 
the thickness of the ring where dust accumulates, estimated as the distance travelled by dust 
particles in the radial direction before being decoupled, SR ~ 3 x 2'7iVr/^, a factor of 3 is 
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assumed for concreteness. The enhancement factor is then readily estimated as 77 ~ (3sZ) ^, 
resulting in the enhanced metallicity in the ring 3e ~ ^73 ~ (SsZq)"^ ^ 1. 

After decoupling, the gas is at a higher temperature and resists gravitational collapse, 
while the dust has smaller velocity dispersion and is thus able to form gravitationally-bound 
objects. The Toomre parameter of dust within the ring decreases inversely proportional 
to dust content in it and, because the dust velocity dispersion Od [pd <^ -Rf^) decreases as 
Gd oc a"^, 

Q. = i^~10A,„(^)\ (9) 

TvGl^d \ a / 

where Oq is the initial radius of the grain and we have assumed that before decoupling dust 
and gas velocity dispersions are equal, ad = Vd- It can be readily shown that during the 
growth stage, the grain radius a increases considerably by factor a/ao ~ 503~^/^m^/'^. This 
means that Qd gradually falls below unity. As a result, dust particles in the ring start 
assembling due to self-gravity, with characteristic free-fall time 



32TTGpd,, 



(10) 



where pd,e is the enhanced dust density in the ring. It can be expressed as pd,e = CSp^Qje, 
where 0.3 .^©36 is the mass fraction of dust in the ring and p ~ mun, and n is assumed to 
be close to the value at the final fragmentation (Omukai et al. 2005), i.e. n ~ 10^3"^. The 
free-fall time is then 

irilS inie 

tffd- J-^ =^ ^^<lQ^h sec. (11) 

It is readily seen that the assembling can take tens to hundreds of rotation periods, depending 
on the mass of the protostellar cloud, rotation free parameter A and metallicity 3. 

It is clear that the overall picture depends on the dust optical model, i.e., on the en- 
richment scenario. When, for example, th e abundance p a ttern of a pair-instability SN with 



a 195 M0 progenitor is applied, for which lOmukai et al.l (j2005[ ) and ISchneider et al.l (|2006[ ) 



predict the density at the final fragmentation step n ~ 10^3"^/^, decoupling condition be- 
comes insensitive to metallicity and dust occupies the whole disk. As a consequence, the dust 
distribution throughout the disk becomes very dilute, making its gravitational accumulation 
a very slow process. At such conditions, however, dust can be selectively assembled by the 
radiation from the central star through the "mock gravity" instability. 
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3.3. Radiation-driven assembling of dust clumps 



When either dust optical properties critically deviate from the model used in above 
estimates, or angular momentum is low A <C 1, the described scenario of centrifugal dust 
accumulation does not operate anymore. In such conditions, a complementary mechanism, 
which selectively drives dust particles to form clumps, connected with radiation pressure 
from a central stellar or protostellar object, can come into play. Let us discuss briefly this 
scenario assuming for the sake of simplicity A ^ 1. When a dusty gas is illuminated by an 
external radiation fleld, a mock gravity ins tability can develop with the growth rate similar 
to the one of Jean's instability (jFieldlll97ll ). 



(12) 



Here Tr = Kc{pr/ pY^"^ where k is the inverse length of photon extinction and c is the light 
speed. Radiation density is pr = inF/c^, where F is the photon energy flux, p = p^ + Pg 
is the total dust and gas mass density, and k is the wavenumber, a = k/n, the normalized 
wavenumber, /(ex) is deflned by 



-a 
3 



1-/3 

(3 



+ 



1 



1 — tan ^ a/a 



(13) 



where (3 is the fraction of incident radiation absorbed by a dust grain. The sound speed 
in dust-gas mixture, Cdg, is essentially equal to the sound speed in gas, c^g — Cg. It can be 
readily seen that for all wavenumbers of interest, a ^1, in which case the asymptotics of 
/(cr) ~ Qa, where Qa is the absorption efficiency (jFieldlll97ll ). In the long-wavelength limit, 
the wavenumbers are less than the critical wavenumber, k < kc = Tr/cdg, and perturbations 
in such dust-gas mixture become unstable against aperiodic growth of clumps under the 
action of radiation pressure. 



The optical flux from a nearly-solar mass protostar at distance r is 







47rr^ 



2 X 10^ r ^ erg cm ^ sec ^ 



(14) 



where r^^ is the distance from star in units of 10 AU and p describes the deviation of the 
luminosity from the main-sequence (MS) solar luminosity before settling of a protostar onto 
the M S. In Eq. f lT^ we explicitly assumed p = 10 following conclusions by iPalla &: Stabler 
(119931 1. We can estimate the inverse extinction length as 



Tia VZlr. 



(15) 
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where dust, as above, is assumed monodispersal. For typical parameters, we obtain the 
critical length 

o o / \ 1/2 o - / \ 1/2 

ATI Zn Cg I P \ Sirpa I p \ Cg 



kc H C \QaPrJ ^^ZqIP \QaPrJ C 

here density of dust particles is assumed as p = 3 g cm~^. 
The corresponding critical mass is 

Here Tioo = T/lOO K is the temperature in the circumstellar cloud and ao.i = a/0.1 pm. 
If one assumes that the density in such a cloud, surrounding the central protostar, is close 
within a factor u < 1 to the density at the time of last fragmentation (when the protostar 
has formed), the critical mass becomes 

~ 1026zy-l/23-9/V g, (18) 

for dust composition related to a ~ IQSM© polluting SN, and 



Mf ~ 3 X lO^V^/^^-Vf^ g, (19) 

for a 22 Mq SN dust composition, as inferred from ISchneider et al.l (120061 ): here Tioo = 0.1 
and p = 10. In both cases, is close to the Jovian mass for rio = 1, when the metallicity 
range 3 = 0.1 — 0.01 is assumed. Note that the corresponding mass of metals, = ZM^ = 
Zq^Mc, is either equal or even less than the mass of the Earth. It is interesting to note 
that the critical mass depends on the distance from the central protostar Mc oc r^, which 
is obviously caused by the fact that the instability is driven by the protostar radiation — 
the farther from the protostar, the weaker is the instability, and the longer is the unstable 
wavelength: at r = 1 AU the critical mass is already in the range of Earth masses. It is 
worthwhile to stress here that low specific angular momentum A ^ 1, that we have assumed 
in this section, does not mean that dust clumps would fall onto the central star and evaporate. 
Firstly, A <^ 1 implies only that the disks form at much later stages than described by Eq. 
([8]), that is, dust does not decouple centrifugally. Secondly, it can be readily seen that the 
characteristic growth time of the mock gravity instability, ~ lOao.irio yr (for a 22 Mq 
dust composition of Schneider et al., 2006), is comparable (or shorter at radii rio > 1) to 
the free-fall time of clumps onto the star, t// ~ lOr^Q^ yr (Eq. ITTl) . When, on the contrary, 
angular momentum is not necessarily A ^ 1, an additional positive energy - epicyclic motion 
(Pringle & King, 2007) - enters the r.h.s. of Eq. f lT2|) along with the sound mode k'^c'^g, 
which increases critical lengths and masses of the fragments formed. 
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4. Conclusions 



1. Recently discovered planets around metal-deficient stars, such as HIP 13044 and HIP 
11952, indicate that planets can form in the early universe. Free-floating planets ob- 
served in MOA experiment can also be old as they are located in the region of stellar 
Pop. 11 (in this sense, they belong to Population 11 planets — planets formed after 
the initial episodes of enrichment by Pop. Ill stars, either in situ with Pop. 11 hosts 
or ex-situ from direct gas collapse). At the same time, it cannot be excluded that a 
well-established increase of metallicity towards the Galactic center can contaminate 
this expectation. 

2. On the other hand, observations deflnitely show that metals in the universe are mixed 
very inefficiently, resulting in a highly inhomogeneous distribution of metallicity in 
the universe: it is clearly seen both in stellar populations of GRB hosts and in the 
IGM clouds of GRB damped Ly-a systems. As a consequence, even in a high redshift 
universe (say, ^ ^ 6), the conditions can exist for formation of planets with normal 
metallicities, either anchored to host stars or free-floating in the star-forming clouds. 

3. Planets can also form in metal-poor environments: 

• In the cases when the speciflc angular momentum in protostellar molecular clouds 
is comparable to the value of the Milky Way molecular clouds. In such case, they 
form preferentially in outer regions of circumstellar disks due to formation of dust 
rings, where dust particles are decoupled from gas and have negligible inward 
radial velocity. 

• If, however, rotation of starforming clouds is small and circumstellar gas accretes 
quasi-spherically, the radiation from the central protostar can stimulate formation 
of clumps through the mock-gravity instability caused by radiation pressure on 
dust grains. The characteristic masses of such clumps for flducial parameters are 
close to the Jovian mass. Such clumps can, in principle, give rise to the formation 
of (free-floating) planets in a metal-deflcient environment of the early universe. 



Therefore, the conditions favourable for planet formation — sufficiently high dust-to-gas ratio 
— could have existed in the early universe in both metal-rich and metal-poor environments, 
which can explain the existence of planets with an age of up to 13-14 Gyr in the local 
universe. These conditions, that we have envisaged here, are however only necessary though 
not sufficient, so that further observational and theoretical study are needed in order to 
firmly substantiate this conclusion. 
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